Intra-Landau-level collective excitations in a bilayer disordered electronic system 
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We investigate intra-Landau-Ievel collective excitations in a bilayer disordered two-dimensional 
electron system exposed to a perpendicular magnetic field. The energy spectrum is calculated 
within the random phase approximation by taking into account electron-impurity scattering in the 
self-consistent Born approximation which includes consistent vertex corrections. Signatures of these 
bilayer excitations in drag and collective excitation measurements are identified. 
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In bilayer two dimensional electron systems (2DES), an 
additional degree of freedom, associated with the layer in- 
dex, produces collective excitations (CE), which have no 
counterpart in the individual 2DES 0|. New interlayer 
electron correlations arise which govern different ferro- 
magnetic and antiferromagnetic phases of the quantum 
Hall states 0j . The bilayer CE also play an essential role 
in the dynamical screening of interlayer Coulomb inter- 
action in the frictional drag effect both in zero H, ^ 111 
and finite magnetic fields H, H S 13 ■ 

The basic source of physical information on the fun- 
damental properties of a system is the spectrum of CE. 
In zero magnetic field, the spectrum of CE in the bilayer 
2DES without disorder effects has been studied compre- 
hensively both in theory and experiment The 
spectrum includes not only the optical wave with in- 
phase density fluctuations and the usual square-root dis- 
persion, but also an additional acoustical branch with 
out-of-phase fluctuations and linear dispersion. In the fi- 
nite magnetic fields, the dispersion of inter-Landau level 
(LL) CE with energies close to the multiples of the cy- 
clotron energy was obtained analytically at the zero tem- 
perature limit [l2] . Experimental measurements of the 
spectrum of inter-LL CE in bilayer 2DES were performed 
only recently: the principal Bernstein magnetoplasmon 
modes were observed by means of inelastic light scatter- 
ing [13. 

Of central importance, however, is the spectrum of 
CE in the extreme quantum limit when no kinetic en- 
ergy change occurs in the 2DES. These CE with ener- 
gies, hu, below the cyclotron energy, Hlub, are referred 
to as the intra-LL CE. In the fractional quantum Hall ef- 
fect (FQHE) regime the intra-LL CE are the well-known 
magnetorotons in single T^l and double layer jlS] 2DES 
with the dispersion minima at wave vectors q ~ 1/£b 
£b is the magnetic length. In disordered 2DES, a 
continuum of particle-hole excitations takes place within 
a broadened LL. In the individual disordered 2DES the 
spectrum of intra-LL CE has been studied by Antoniou 
et al. 01 for the lowest LL and it has been found that 
the plasmon-like CE emerge above the particle-hole con- 
tinuum (PHC). It has been also predicted by Merkt 
that in the disordered 2DES, electron-electron (e-e) in- 



teraction overcomes Kohn's theorem [l9| and results in 
a low-lying disorder-induced cyclotron resonance mode. 
Recently, the importance of intra-LL CE in the drag ef- 
fect has been also discussed 

In the present paper we calculate the spectrum of intra- 
LL CE in the bilayer disordered 2DES in a perpendicular 
magnetic field at finite temperatures when a higher LL is 
partially filled. We treat e-e interaction within the ran- 
dom phase approximation (RPA) and electron-impurity 
(e-i) scattering in the self-consistent Born approximation 
(SCBA) [23, which includes the consistent vertex cor- 
rections to the response function of the individual 2DES 
[itI . We obtain that the whole spectrum of intra-LL CE 
is located in a finite wave vector interval and, depending 
on the system parameters, there can exist forbidden 
zones within this wave vector interval (Iq is the Landau 
index of the outermost partially flUed LL). Although our 
approach allows us to take into account coupling of the 
broadened LL at finite temperatures, the coupling is weak 
in the regime of our interest here. Therefore, the response 
function becomes small and the 2DES unable to respond 
in the forbidden wave vector zones around qr where qr 
are determined from the zeros of the Laguerre polyno- 
mial, Li^{tr) = and U = {qrlsf 1^ ir = l,2,...,Zo)- In 
each allowed wave vector zone, we obtain the spectrum 
of bilayer intra-LL CE which consists of two excitations 
with in-phase and out-of-phase density fluctuations and 
with the energies depending on the interlayer spacing and 
the characteristics of individual 2DES. 

The bilayer dielectric tensor (g, uj) (i, j = 1,2 are the 
layer indices) is given in terms of the irreducible electron 
polarization function Hij(g, w), which is obtained from 
the solution of a matrix Dyson equation for the dynam- 
ically screened Coulomb interaction V{q,uj). It has the 
following usual form 

Viq,co)=viq) [l~viq)fl{q,Lo)Y' =v{q)e-\q,co) (1) 

where all quantities are 2x2 matrices, v denotes the 
bare Coulomb interaction tensor. For simplicity we as- 
sume that the layers are infinitely thin and use vn = 
V22 — V ~ 27re^/Kog and V12 = V12 = cxp(— gA)!;, 
where kq is the static dielectric constant and A the in- 
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terlayer spacing. We neglect tunneling between the lay- 
ers and adopt an approximation where the interlayer 
polarization functions are zero, II12 = TI21 = 0. Then, 
the components of the inverse dielectric tensor are given 
as erM2(9,w) = e2,i(q,w)/£bi(g,w) and e^^^^^{q,uj) = 
V12 iq,u)) /shi (<7, w). Here, £1,2 (<7,c^) = l-uiinii^22('?, t^) 
and nii^22('?,^) are the screening and polarization func- 



tions in each layer. The bilayer screening function 
ebi('Z,i^) is the determinant of the dielectric tensor, 
ebi('?,'^) = ei(g,w)e2(g,w) - t;i2(g,w)nii(g, w)n22(g,t^). 

The RPA approximates the exact irreducible II{q,Lu) 
by its simple bubble diagram with respect to the e-e in- 
teraction, which, for finite magnetic fields, gives 



1 f dE { 

n(9,^) = ^ E W / —fT{E - EF)lTn GfiE) j- 
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G^iE + hw) 



[1 - 5-f^,"it, E,E + f^)] [1 - E,E + huj) 



G^iE-huj) 
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(57,f (t, E,E±hio)^ -fQw (t) G1{E)G\, {E ± Hlo). 



For brevity we suppress the layer indices here, 
G" ' 



R,A, 



J {E) are the electron retarded and advanced 

Green's functions, dressed by e-i interaction. In the 

SCBA they are diagonal and given by Gf'^' 

1 



2{E-Ei + d{E-Ei) 



J- 



\E) ^ 
with the imaginary 



part of the square-root taken positive (negative) for the 
advanced (retarded) function. In the short-range im- 
purity model, the half-width Fq — \J ^Hlub^ of the 
El — (I + l/2)hujB LL is independent of the Landau 
index / (r is the transport relaxation time, determined 
from mobility /i) . The bare vertex functions Qui (t) 
are given by the gauge invariant part of the in-plane 



form factor: Qw (t) ^ (-l)'+''e^*i, (t) L\7^' (t) where 
t = [qlsY /2 and L\ [t) is the associated Laguerre poly- 
nomial. The denominators in the r.h.s. of Eq. (j2J differ 
from one due to the vertex corrections 5^f^ which are 
obtained consistently in the short-range impurity model 
with the assumption that the LL are clearly resolved . 
The chemical potential Ep{n, B, T, n) in the Fermi distri- 
bution function fx is determined implicitly by the elec- 
tron density via n = — {2Tr^£g) ^ S;^o /o°° dEfriE — 
EF)l'aiGi{E). Here we use the electron Green's func- 
tions which correspond to the Gaussian density of states, 
lmGi{E) = V2^/To exp (-2 (E - Eif /F§) , without 

unphysical edges of the Landau band We assume 

spin degeneracy so the capacity of each LL is doubled. 

It is seen in Fig. ^ that Re e{q,u)) of the individ- 
ual 2DES becomes negative for intermediate values of 
Lo of the order of Fq. This requires two crossings of 
the Re e(g,a>) ~ axis. The low-frequency zeros 
(LFZ) of Re e{q,u}) always lie within the PHC where 
Im £(q,Ld) is large. The high-frequency zeros (HFZ) are 
close to or above the upper edge of the PHC where 
lTCLe{q,Lo) is small or exactly zero within the SCBA. 
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FIG. 1: The structure factor (solid line), the real (dashed 
line) and imaginary (dotted line) parts of dielectric function 
of the bilayer 2DES with the A = 30 nm spacing and for 
lo — 4. The thin lines correspond to the individual 2DES in 
which Im e{q,uj) differs from that of in the bilayer 2DES by 
a numerical factor and is not shown here. 



The pronounced peaks of the dynamical structure factor 
S{q,uj) oc — Im (l/e(q, oj)) as a function of lu correspond 
to HFZ, while no peaks correspond to LFZ. Therefore, 
HFZ represent intra-LL CE of the system while LFZ can- 
not be interpreted as excitations. We believe, however, 
that these particle-hole states play an important role, 
particularly in realization of the frictional magnetodrag 
effect, and we calculate LFZ and present them below in 
the spectrum together with the intra-LL CE modes. 

In the bilayer 2DES, there appears a pair of zeros in- 
stead of LFZ and HFZ of the single layer 2DES and 
Re e{q,u!) becomes negative only within small regions, 
restricted by the two zeros in each pair. Therefore, the 
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FIG. 2: The dispersion relations of intra-LL CE referred to 
the lo = 4: half-filled LL in the individual 2DES for the three 
values of mobility. 



FIG. 3: The dispersion relations of intra-LL CE in the first 
allowed wave vector zone for the bilayer 2DES. The thin lines 
represent the spectra in the single-layer limit. 



bilayer intra-LL CE spectrum consists of two branches 
which correspond to the in-phase and out-of-phase intra- 
LL CE. Accordingly, one can see from Fig. ^ that the bi- 
layer structure factor, instead of the single peak of S{q,Lj) 
of the individual 2DES, shows a double-peak structure 
near the upper edge of the PHC, while again it demon- 
strates monotonic behavior at low energies. 

In Fig. 121 we plot the spectra of intra-LL CE of individ- 
ual 2DES at finite temperatures and for various values of 
electron mobility. Inclusion of the consistent vertex cor- 
rections allows us to reproduce the square-root like dis- 
persion ^3 of intra-LL CE at small momenta. However, 
no exact q limit exists in the spectrum. In the dis- 
ordered 2DES the density of states remains finite in the 
centrum of the LL and the response function vanishes as 
q^, even if oj is in the immediate vicinity oi ujb- In the 
short wavelength limit, again the disordered 2DES does 
not respond since the response function vanishes expo- 
nentially at wavelengths much shorter than £b- When 
the lowest LL is partially filled, there exists only one 
allowed wave vector zone, and our findings for = 
are in agreement with Ref. For CE referred to the 
higher LL, one can see from Fig. El that the number of 
allowed wave vector zones can be multiple and increases 
with II. Near half-filling of the = 4 LL and at T = 0.1 
K there exist only one, two, and three allowed zones for 
/i = 35,200, and 1000 T^^, respectively. Notice that for 
Iq = A the first two of the maximum four available wave 
vector gaps are located around g^^B = 0.8 and 1.87. In 
the first zone close to g = 0, for all three values of fi, the 
intra-LL CE emerge above the PHC in the middle of the 
zone where the energy shows a maximum. In this region 
the spectrum describes CE modes without damping and 
delta-function like peaks occur in S{q,uj). Towards the 
edges of the allowed zone the energy decreases. Within 



the PHC hu! < 2To near its upper edge, the peaks of 
S{q, uj) become broadened, and though the intra-LL CE 
modes are damped but still are well defined excitations. 
Exactly at the edges of the allowed zone, the group ve- 
locity shows anomalous behavior and becomes infinite. 
At these points CE merge with the particle-hole states, 
which are described by LFZ of Re s{q,ui). Below this 
merging energy the damping is so strong that the zeros 
of Re e{q,u)) cannot be interpreted as CE. The greatest 
maximum of energy and the largest width of the allowed 
wave vector zone are achieved, when the outermost LL 
is near half-filling, /i large, and T low. In the higher 
zones the energy maximum is smaller and the intra-LL 
CE emerge above the PHC only in the second zone for 
the highest value of mobility. 

In Fig. 13 we plot the dispersion relations of intra-LL 
CE of the bilayer disordered 2DES with the interlayer 
separation A = 30 nm. We restrict our consideration 
here to the symmetric bilayer 2DES with matched elec- 
tron densities and present calculations only for the first 
allowed wave vector zone. In this case £1 = £2 = £ and 
the bilayer screening function ehi{q,io) is represented as 

£bi(g, cj) = (1 - a^) {e{q, uj) - £in(g)) {£{q, t^) " £out('?)) 

(3) 

where we have defined fin.outCg) = 1/(1 ± a ^) with 
a — exp(— gA). Now the spectrum of intra-LL CE in the 
bilayer 2DES is determined implicitly by the conditions 
Re £(g, Win^out (g)) = ein,out(g)- 

As seen from Fig. |3| the spectrum of bilayer intra-LL 
CE is split into the in- and out-of-phase modes. The 
in-phase CE mode, irrespective of the spacing A, is less 
affected by the interlayer interaction, lies above and be- 
haves similarly to the CE mode of individual 2DES. The 
out-of-phase CE branch lies below the single layer mode 
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and has a dipole oscillator strength. At small wave vec- 
tors the dispersion of in-phase modes is close to the 
square-root one while the out-of-phase modes demon- 
strate a smoother dispersion and achieve the maximum 
energy at larger values of q. The splitting energy, SE, 
decreases in A which supports the tendency that within 
the limit of A —^ oo, the dispersion curves of in- and 
out-of-phase modes degenerate and coincide with that of 
single layer intra-LL CE mode. SE has a minimum near 
the upper edge of the PHC and close to the right edge 
of the allowed wave vector zone. Towards the right edge 
of the zone, 6E increases slightly, while the increase is 
strong towards the left edge. The CE modes with q close 
to the edges of the zone and with energies huj < 2To are 
responsible for the enhancement of the frictional magne- 
todrag effect. The out-of-phase CE mode has relatively 
strong damping and contributes largely to the drag ef- 
fect. However, in certain ranges of wave vectors close to 
the edges of the allowed zone, the in-phase modes develop 
and exist solely and will contribute to the drag effect in 
the absence of out-of-phase modes. 

The corresponding bilayer structure factor for the Iq = 
4 LL is depicted in Fig. ^ which we calculate in the 
plasmon-pole approximation. For qls = 0.33 the split- 
ting energy is approximately equal to O.lTFo. For /i = 35 
T-i and B « 3.22 T we have To ~ 4.84 K and this gives 
SE more than 0.8 K. This is a quite measurable quan- 
tity, however, notice that at smaller wave vectors, SE is 
appreciably larger (Fig.O. As seen from Fig.Q] the out- 
of-phase peak of S(q, uj) is about six times lower than 
the in-phase mode peak. However, both the in- and out- 
of-phase CE modes have enough weight to be observed 
experimentally, and the signatures of these bilayer intra- 
LL CE in the infrared absorption and inelastic light scat- 
tering measurements should constitute asymmetric dou- 
blets. 

We expect that the approximations we use here are 
adequate for the 2DES in which disorder dominates over 
Coulomb correlations. This is safer in the first allowed 
zone where g/s < 1 and the RPA is an acceptable ap- 
proximation even in the disorder free case \m. In the 
disordered 2DES, the RPA is vahd when no FQHE exci- 
tation gap is developed 17] . This is the case if the gap in 
the disorder free case, Aq, is smaller than the typical dis- 
order energy According to the Laughlin's estimates 
m Ao « 0.056Ec for the smallest fraction v = 1/3, 
for which Aq is maximal. Here, we consider moderate 
magnetic fields well below the FQHE threshold |2^ and 
obtain the intra-LL CE referred to the higher LL. For 
/J, = 35 T^^, the ratio Ec/Tq is less than 20, and near 
half-filling of the Iq = 4 LL, the system is well away from 
the regime where FQHE occurs. Thus, we expect that 
the regime considered here is favorable for the applicabil- 
ity of this theory and for the observation of these novel 
bilayer intra-LL CE with attractive many-body and dis- 
order effects. 
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